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Abstract 

We treat analytically a model that captures several features of the phenomenon of spatially 
inhomogeneous reversal of an order parameter. The model is a classical Ginzburg-Landau 
field theory restricted to a bounded one-dimensional spatial domain, perturbed by weak spa- 
tiotemporal noise having a flat power spectrum in time and space. Our analysis extends the 
Kramers theory of noise-induced transitions to the case when the system acted on by the noise 
has nonzero spatial extent, and the noise itself is spatially dependent. By extending the Langer- 
Coleman theory of the noise-induced decay of a metastable state, we determine the dependence 
of the activation barrier and the Kramers reversal rate prefactor on the size of the spatial do- 
main. As this is increased from zero and passes through a certain critical value, a transition 
between activation regimes occurs, at which the rate prefactor diverges. Beyond the transition, 
reversal preferentially takes place in a spatially inhomogeneous rather than in a homogeneous 
way. Transitions of this sort were not discovered by Langer or Coleman, since they treated 
only the infinite-volume limit. Our analysis uses higher transcendental functions to handle the 
case of finite volume. Similar transitions between activation regimes should occur in other 
models of metastable systems with nonzero spatial extent, perturbed by weak noise, as the size 
of the spatial domain is varied. 

1 Introduction 

The phenomenon of noise-induced escape from a metastable state, and the related phenomenon of 
noise-induced transitions between the stable states of a bistable system, occur in many places in 
the sciences and engineering [ 1 1. Our focus in this paper is the case when the system on which the 
noise acts has nontrivial spatial extent, and the noise is not spatially uniform, but has significant 
spatial dependence. This case arises in numerous physical contexts, in both equilibrium and non- 
equilibrium systems. Examples of noise-induced transitions in spatially extended systems include 
pattern formation in convective and electroconvective systems far from equilibrium |El|3l, ther- 
mally activated magnetization reversal in nanomagnets 0J|3, vortex creep in superconductors 0, 
the growth of instabilities in metallic nanowires \7\, and many others. 
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Possibly the simplest mathematical model of a spatially extended system undergoing a noise- 
induced transition is a Ginzburg-Landau scalar field theory perturbed by weak spatiotemporal 
noise. Even here, the mathematical complexity of the problem has restricted most treatments to 
the case of one spatial dimension. The phenomenon of classical nucleation on a line was treated by 
Langer in a seminal paper, and a closely related quantum phenomenon was treated by Callan 
and Coleman [ 9 1 . In the limit as the size of the one-dimensional spatial domain tends to infinity, 
they were able to compute (respectively) the nucleation rate per unit length, and what amounts to 
a tunnelling rate per unit length. Their work has been widely applied and cited. 

However, many practical applications require an understanding of the phenomenon of noise- 
induced nucleation in systems with a spatial extent that is both nonzero and non-infinite, and 
the extent to which the nucleation rate per unit length depends on system size. Noise-induced 
transitions between the two stable states of a bistable nonlinear field theory, in a finite-volume 
domain, have been investigated by several authors ifTUl ITT1 IT21. The difference between finite 
and infinite systems is not merely quantitative. Recently, the authors 1~i~3l uncovered an unusual 
effect, akin to a phase transition, that occurs in an overdamped classical Ginzburg-Landau field 
theory with a bistable <p A potential, as the length L of its one-dimensional spatial domain is varied. 
A similar phenomenon had previously been seen in a quantum context ifbH [T51 . The lowest- 
energy 'saddle' between the two stable configurations, through which noise-activated transitions 
preferentially occur (especially when the noise strength is low) may bifurcate. Below a critical 
length L c , this transition state is a spatially constant field configuration, but at L = L c it bifurcates 
into a spatially varying pair of configurations, degenerate in energy. Subsequent work ifToTl shows 
that a similar bifurcation occurs in a Ginzburg-Landau model with an asymmetric potential, in 
which the two new spatially varying transition states are not degenerate; so one or the other is 
preferred. The 'phase transition' at a critical length L c is therefore reasonably robust. 

So at a critical system size, there may occur a major change in the phenomenology of noise- 
induced transitions between the stable states of a Ginzburg-Landau field theory, which proceed 
via nucleation. Associated with the bifurcation of the transition state is a bifurcation of the MPEP 
(most probable escape path, i.e., transition path in configuration space) extending uphill from each 
stable state to the transition state. The bifurcation is driven by this preferred nucleation pathway 
becoming unstable in the transverse direction. As one would expect, the transition rate is strongly 
affected by the bifurcation. Formally, the prefactor in the Kramers (weak-noise) nucleation rate 
diverges at L = L c . This signals that precisely at L = L c , the phenomenon of transition between 
the two stable states becomes non-Arrhenius: the rate at which it occurs falls off in the limit of 
weak noise not like an exponential (with a constant prefactor), but rather like an exponential with 
a power-law prefactor. 

Our paper Ref. lfT3ll computed the Kramers prefactor as a function of L, through the bifurcation, 
by using elliptic functions. The use of higher transcendental functions made it possible to treat the 
case L < oo, i.e., to go beyond the analyses of Langer, and Callan and Coleman. However, we 
considered in detail only the case when Dirichlet conditions are imposed on the Ginzburg-Landau 
field at the endpoints of the spatial domain. How sensitive is the occurrence of a phase transition at 
some L = L c to the boundary conditions employed? We examine this question here, by extending 
our qualitative results and nucleation rate computations to the cases of Neumann, periodic, and 
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Figure 1 : Bistable potential for the order parameter given by Eq. [TJ 

antiperiodic boundary conditions. We shall show that the phase transition is robust to changes in 
boundary condition, but that unsurprisingly, there are quantitative differences among the various 
cases. Examination of these differences will permit us to present our methods in more detail, 
and provide further insights into the behavior of noise-activated transitions in spatially extended 
systems. 



2 Model and Phenomenology 

Consider an order parameter governed by a bistable quartic potential energy function that may 
be written in dimensionless units as 

V(0) = (0 2 -l) 2 /4. (1) 

This potential is shown in Fig. [lj The stable states = ±1 are degenerate, with zero energy. If 
evolves in a deterministic, overdamped way, it will satisfy the equation = — V '((/)). If it is 
perturbed by additive white noise, it will satisfy, instead, = — V (0) + e 1 / 2 ^(t). Here £ is unit- 
strength temporal white noise, satisfying {£{ti)£(t 2 )) = 5(ti — t 2 ), and e is the noise strength. In 
a thermal context, e oc kT. 

We shall treat here the case when is actually a classical field on a one-dimensional spa- 
tial domain: in dimensionless units, the interval [0, L]. The simplest extension of the overdamped 
stochastic evolution equation to incorporate spatial effects is the stochastic Ginzburg-Landau equa- 
tion 

= cj>" - V(4>) + e 1 / 2 ^) (2) 

= 0" + 0-0 3 + e V2 e(X;t)) 

where £(x, t) is unit-strength spatiotemporal white noise, which is defined to satisfy (£ {x\, ti)£ (x 2 , t 2 )) = 
5(xi — x 2 )5(ti — t 2 ). Zero-noise dynamics will be 'gradient', i.e., conservative. That is, if e = 0, 
then 

= -mm , O) 
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where 
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is the energy functional. So the statistical properties of the stochastically evolving field <p are de- 
scribed by equilibrium statistical mechanics. However, just as in the model without spatial extent, 
nonzero noise can induce transitions between the stable states, which naively are the two field con- 
figurations 4> = ±1. This will typically occur via nucleation. A droplet of one stable configuration 
will form in a background of the other, and will 'nucleate' : under the driving influence of the noise, 
it will spread to fill the entire spatial domain. Of course, it is more likely for a small droplet to 
shrink and vanish, especially when the noise is weak. 

In the infinite-dimensional configuration space, an MPEP goes 'uphill' from each stable field 
configuration, leading to a preferred transition configuration (saddle) that lies between them. Each 
of these MPEPs is a nucleation pathway, i.e., a path of least resistance. By time-reversal invariance, 
each MPEP is a time-reversed zero-noise 'downhill' trajectory. If the noise is weak, the order 
parameter is expected to flip between the two stable configurations in a Markov way, with the 
expected waiting time in the basin of attraction of each being an exponential random variable, as is 
typical of slow rate processes. The activation rate (the reciprocal of the mean time between flips) 
will be given in the e — > limit by the Kramers formula 



Here AW is the activation barrier, which quantifies the extent to which the preferred transition 
configuration between the two stable configurations is energetically disfavored, and T is the rate 
prefactor. Due to the normalization convention implicit in ©, AW = 2AE, where AE is the 
energy of the transition state minus the energy of either stable state. Here energy is computed 
from the functional 7i[-] of ©. The factor of 2 arises from our decision to multiply unit-strength 
spatiotemporal noise by e 1 / 2 rather than by 

The calculation of AW and T will make up the technical component of this paper. Since we 
are using dimensionless units, they will depend only on the length L and the choice of boundary 
conditions at the endpoints x = and x = L. The boundary conditions affect the way in which 
order parameter reversal occurs, since they may force nucleation to begin, preferentially, at the 
endpoints. It should also be noted that for some choices of boundary conditions, the stable states 
may only be approximations to the uniform </> = ±1 configurations. 



We begin by examining the stationary solutions of the noiseless (e = 0) evolution equation. These 
will depend on the interval length L and the boundary conditions. The four boundary conditions we 
consider are periodic (P), antiperiodic (AP), Dirichlet (D), and Neumann (N). All four, applied to 
the stochastic Ginzburg-Landau equation, have potential applicability in physical modelling. For 
example, in modelling thermally activated magnetization reversal in a finite-length ferromagnetic 
nanowire, Neumann boundary conditions are appropriate ifTTI . The original treatments of Langer 
and Callan-Coleman used periodic boundary conditions; and so forth. 



r ~r exp(-Aiy/e). 



(5) 



3 The Stable and Transition States 
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Periodic. — The conditions are 0(0) = 0(L), 0'(O) = 4>\L). There are three constant time- 
independent solutions: = il> eacn with energy 0; and = 0, with energy 1/4. It is easy 
to see that = ±1 are stable for any L, and = is always unstable. However, nonconstant 
solutions may also exist for a range of L. These will be considered below. 

Antiperiodic. — The conditions are 0(0) = — 0(L), 0'(O) = — 0'(L). The only constant time- 
independent solution is = 0. 

Dirichlet. — Only the conditions 0(0) = 4>(L) = will be considered here. Many other Dirichlet 
boundary conditions could be examined, but this case is particularly interesting because when L is 
sufficiently large, it has two nonconstant stationary solutions, which permits the possibility of a 
noise-induced transition. This is the case that was treated in Ref. lfi"3ll . 

Neumann. — Though there are many possibilities, we consider here only the conditions 0'(O) = 
0'(L) = 0. The constant solutions here are the same as for periodic boundary conditions, but it 
will be seen that there are important differences between the two cases. 

What nonconstant time-independent solutions exist when e = 0? By ©, any stationary solu- 
tion satisfies 

0" = -0 + 3 , (6) 

subject to the given boundary condition. That is, = —5H/5(j) = 0. The linearized noiseless 
dynamics in the vicinity of such a state are specified by the Hessian operator 5 2 H/5(f) 2 . Sta- 
tionary configurations for which this operator has all positive eigenvalues are stable; those with 
a single negative eigenvalue are possible transition states. For the latter, the eigenvector corre- 
sponding to the negative eigenvalue is the direction along which the MPEP approaches, in the 
infinite-dimensional configuration space. Nonconstant transition states are often called instanton 
states, in a nomenclature derived from Callan and Coleman S- In 4 and other theories, the term 
'instanton' usually refers to a kink-like field configuration = 0(x) asymptotic to = ±1 as 
x — > ±oo, with a single node lfi"8l . 

To satisfy the specified boundary conditions on [0, L], a rather different sort of instanton state 
must be used. It is easy to check that the so-called periodic instanton = 4>- ms t, m (x) is a time- 
independent solution of © for any m in the range < m < 1. Here (cf. Ref. EDX 



0inst,m(^) = \ — - sn(x/ V m + 1 I m) , (7) 

V m + 1 

where sn(- | m) is the Jacobi elliptic sn function with parameter m. Its quarter-period is given 
by K(m), the complete elliptic integral of the first kind [20], which is a monotonically increasing 
function of m. It can be viewed as an infinite alternating sequence of kinks and anti -kinks, spaced 
a distance 2K(m) apart. As m — > + , K(m) decreases to n/2, and sn(- | m) degenerates to 
sin(-). As m — > 1~, the quarter-period increases to infinity (with a logarithmic divergence), and 
sn(- | m) degenerates to the nonperiodic function tanh(-), which is the canonical single-kink 
sigmoidal function. It is no accident that the hyperbolic tangent function appeared in the Langer 
and Callan-Coleman analyses in connection with the limiting (L — > oo) shape of a critical droplet. 
For a careful recapitulation, see Ref. I 2"T1 . 

In the present context, the value of m in © is determined by the interval length L and the 
boundary conditions, as follows. 
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Periodic. — As stated earlier, the uniform configurations = ±1 are the stable states for all L. If the 
nonuniform solution © is to satisfy the boundary conditions, then L must be an integer multiple 
of a full period. It is obvious that larger integers correspond to higher energies (and therefore 
activation barriers); so the physically relevant nonuniform transition state is = (j>xDst,m uP {x), 
where m U}P is determined implicitly by 



The subscript u denotes an unstable (i.e., transition) state and P denotes periodic boundary condi- 
tions. Similarly, s will denote a stable state, and the other boundary conditions will be correspond- 
ingly abbreviated. 

This nonuniform stationary state is not present for all L: there is a minimum value of L, 
corresponding to m = 0, below which no solution of the form = 0i ns t,m u P {x) can be constructed. 
That is, because K(0) = n/2, the nonuniform transition state cannot occur below L = 2n. Below 
this value, the only possible transition state is the uniform state = 0. This solution remains an 
unstable stationary state for L > 2ir, but with higher energy than that of the instanton state. 

These considerations suggest the following qualitative picture of the periodic-b.c. case, to be 
justified subsequently by an energy and eigenvalue analysis. Suppose the system is in one of 
the two stable states = ±1. How does a transition induced by (weak) noise to the other state 
proceed? If L < 2n, the system preferentially 'climbs uphill' to the = state, after which it 
'rolls downhill' to the other stable state. So reversal of the order parameter preferentially takes 
place in a spatially homogeneous way. If L > 2n, however, there is another saddle, lower in 
energy than the instanton state and therefore a much more probable intermediate state when the 
noise strength is small. This is the instanton configuration = 0i ns t,m u P { x )i which is positive on 
half the interval and negative on the other half. This configuration can be thought of as a 'droplet 
pair', in which half the spatial domain is occupied, to a first approximation, by each of the two 
stable values of 0. Equivalently, it can be viewed as a configuration with a kink and an anti-kink, 
or two Bloch walls (if one is modelling a magnetic system). A transition occurs in the following 
way: a small droplet of the reversed value for the order parameter is formed, and under the driving 
influence of the noise, it grows to occupy occupy half the interval. At that point, the saddle has 
been reached, and the droplet continues growing of its own accord until it fills the remainder of 
the interval, i.e., 'rolls downhill'. Both halves of the transition can be interpreted in terms of kink 
motion: the uphill half being noise-driven, and the downhill half being deterministic. 

Actually, there is a complication here: due to the translation symmetry that accompanies 
periodic boundary conditions, the transition state is necessarily infinitely degenerate. That is, 
= 0inst,m u P (x — x ), for any x , will serve as a transition state. In physical language, the kink- 
anti-kink pair, i.e., the droplet, may form anywhere along the interval, leading (in the L — > oo 
limit) to a transition rate per unit length. We treat this matter elsewhere 1221 (see also Ref. ifToTl ). 
Antiperiodic. — In this case there is a nonuniform stable state of the form = 0mst,m sAP , with 
m S) AP defined implicitly by 



However, this equation has no solution when L < it; in this range there is a single stable state 
given by = 0, and no possibility of an accompanying transition. The = solution becomes 




(8) 




(9) 
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Figure 2: The stable states (S) and transition states (U), for periodic (P), Dirichlet (D), and Neu- 
mann (N) boundary conditions, when L — 10. Each nonzero state has a degenerate counterpart 
obtained by <^> i — ^ — 0, and UP may be shifted arbitrarily, so it is infinitely degenerate. (From 
Ref. El.) 



unstable when L > n, and in fact is the transition state for n < L < 3n. For larger L this transition 
state undergoes a pitchfork bifurcation into a pair of instanton states with three half-wavelengths 
in the interval [0, L). That is, 

6\Jm UiAP + 1 K(m UyAP ) = L . (10) 

As in the periodic case, this nonuniform transition state is infinitely degenerate; and the nonuniform 
stable state is too. 

Dirichlet. — This case was examined in detail in Ref. 1131 . If L < ir. the model is monostable with 
the configuration = as the only stable state. When L increases through tt, this state undergoes 
a pitchfork bifurcation into a pair of stable configurations with 

2^ + lKK D ) = L. (11) 

In the range < L < 2n, the = configuration is the transition state. It undergoes a pitchfork 
bifurcation at L = 2n; for larger values of L the parameter m of the transition state is determined 

by 

^m u ,D + 1 K(m UjD ) = L . (12) 

Note that the Dirichlet condition (and the Neumann condition, to be discussed below), pin the 
nonuniform stationary states at the boundaries, so that there is no issue of infinite degeneracy 
arising from translation symmetry to complicate the analysis. In physical terms, each of the two 
preferred nucleation pathways begins with the formation of a droplet at either end of the interval. 
As L — » oo, it is not natural to speak of a transition rate per unit length. 
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Neumann. — Here the stable states are the uniform states <f> = ±1. Also, if L < ir the transition 
state is the uniform configuration = 0. When L is increased through ir, the transition state 
bifurcates into a pair of instanton states with the first argument of the Jacobi elliptic sn function 
in © shifted by K(m), a quarter wavelength (equivalent here to L/2). That is, the transition state 
is given by 

sn(x/Vm + 1 + K(m) | m) , (13) 



V m + 1 

with m = m U:N , a quantity determined implicitly by 



2^fm U)N + 1 K(m u>N ) = L . (14) 

So m U: N = m s ,D- In fact, up to a uniform shift, the Neumann-case transition state is the same as 
the Dirichlet-case stable state. 

The results so far are summarized in Fig.|2j which depicts the stable and transition states for the 
cases of periodic, Dirichlet, and Neumann boundary conditions (the antiperiodic case is excluded 
for figure clarity). Each plotted nonconstant curve consists of an integer number of kinks or anti- 
kinks; or more accurately, an even number of half-kinks or half-anti-kinks. 



4 The Activation Barrier 



As mentioned in Sec. |3J the exponential falloff of the transition rate in the limit of weak noise, 
i.e., its Arrhenius behavior, is determined by the activation barrier AW between the two stable 
states. We noted that AW = 2AE = 2(H[(f) u } - H[(p s ]), with the energy functional H[4>] given 
by @. We shall call AE the 'activation energy'. The calculation of H.[(j)J\ and H[(p u ], the stable 
and transition state energies, is trivial in the case of uniform states, and reasonably straightforward 
in the case of nonuniform (instanton) states, if standard elliptic function formulas |20| are used. 
The results are as follows. 

Periodic. — If L < Lf = 2n, then trivially, AE = L/4. If L > Lf , then 



AE 



1 



3^1 + "Vp 



8E(m u p) - 



(1 - m Ut p)(3m u , P + 5) 



K(r?v 



(15) 



(1 + m UtP ) 

where E(m) is the complete elliptic integral of the second kind ll20ll . Note that as m u iP — > + , 
i.e., L — > (27r) + , the two regions connect in a smooth manner: AE is is differentiable with respect 
to L at L = 2ir, though not twice differentiable (the same relatively smooth join will obtain in all 
succeeding cases). 

As m U) p — > 1~ (i.e., L — > oo), AE — > 4a/2/3. This is a familiar quantity: it is twice the energy 
of a 4 kink lfi"8ll . The factor of 2 is due to the transition state configuration U , the stationary 
periodic instanton solution, having two nodes, i.e., two swings between = ±1 as x varies from 
to L. In field theory language, it comprises a kink and an anti-kink. 
Antiperiodic. — In the range n < L < 3n, 

2 



AE 



3(l + m s , AP )3/2 



[m S)AP + 2)K(m sAP ) - 2{m sAP + l)E(m s , AP ) 



(16) 
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Figure 3: The activation energy AE as a function of the interval length L, for periodic (P), 
antiperiodic (AP), Dirichlet (D), and Neumann (N) boundary conditions. Bullets indicate critical 
interval lengths, at which bifurcations take place (L^ = ir; = L P = 2ir; and Lf p = 3ti). 



If L > L^ p = 3vr, then 
3 



AE 



yfl + {m uAP ) 

1 

~ 6 



Efm, 



1 



,AP 



1 + (m s Ap) 



=E(m„ 



AP 



(17) 



3(1 - m uA p) 
(l + m u , AP )V* 



m S Ap 



;i + m^ P )3/2 



(3m S:AP + 5)K(m SjAP ) 



It is noteworthy that as m Uj AP —> 1 _ , AE — > 4\/2/3, just as in the periodic case. In the periodic 
case, the transition state has two nodes while the stable state has zero nodes; in the antiperiodic 
case, the transition state has three nodes while the stable state has one node. It is therefore not 
surprising that in both cases, the energy difference converges to twice the energy of a <f) 4 kink 
as L — > oo. Note that for any finite L, the relevant m values differ between the periodic and 
antiperiodic cases, but all approach unity from below as L — > oo. 

Dirichlet. — If L < Lf = 2-k, then 



AE 



3(l + m s , D )3/2 
which has the same form as (fl~6b . If L > Lf , then 

2 1 



m SjD + 2)K(m S)£) ) - 2(m S)D + l)E(m S)D ) 



(18) 



AE 



=E(m, 



1 + (m u . 



D 



1 + (m s n) 



=E(m 



(19) 



2(1 - m UtD ) 

(l+m u n)3/2 



(3ra u ,D + 5)K(m UjD ) 



1 - 

(1 + m s , D )3/2 



(3m s , D + 5)K(m S)D ) 
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In the L — > oo limit, AE — > 2\/2/3, which is half the value for the periodic and antiperiodic 
cases. This is because the transition state now has only one more node than the stable state, rather 
than two. The same is true of the Neumann case. Hence these two latter cases will have the same 
activation energy in the L — > oo limit, just as the periodic and antiperiodic values for AE converge 
in this limit; and it will be the energy of a single </> 4 kink. 

Neumann. — If L < = it, then AE = L/4; both the stable and transition states are the same 
here as in the low-L regime of the periodic case. If L > L^, then 



AE 1 



4(1 + m u , N )E(m UtN ) - -(1 - m UiA r)(3m„ iiV + 5)K(m U)P ) 



(20) 



3(l + m Mi7V )3/2 

As noted, the activation energy in the L — ► oo limit equals 2 a/2/3, as in the Dirichlet case. 

The above formulas for the activation energy AE as a function of L are plotted in Fig. [3] The 
bifurcations at L = ir, L = 2n, and L = 3ir are apparent, as is the differentiability (and lack of 
twice differentiability) through each bifurcation. 



5 The Transition Rate Prefactor 

As explained, the activation energy is (up to a factor of 2) equal to the activation barrier AW 
in the Kramers transition rate formula T ~ r exp(— AW/e), e — > 0. Calculation of the pref- 
actor r is a much more involved matter. In general, it requires an analysis of the transverse 
fluctuations about the MPEPs that go uphill to the transition state. In any 'zero-dimensional' (i.e., 
non- spatially-extended) equilibrium system, in which the noise-perturbed evolution equation is a 
stochastic ordinary differential equation, it is known that it suffices to compute the eigenvalues of 
the linearized dynamics at the endpoints of each MPEP, i.e., in the vicinity of the stable states and 
transition state ll2"3l . For example, a system with a two-dimensional order parameter <p satisfying 
the overdamped evolution equation = u(0) + e 1 / 2 ^, where u = u(0) is a drift field that is the 
negative gradient of a potential, will have the Kramers escape rate 



r ~ 2^I A| H A||(S) \ TJV) ex P(- AW V £ ) , e-0. (21) 

Here S denotes the stable fixed point of u and U the saddle point over which escape from the 
basin of attraction of S preferentially occurs. X\\(S) = —dus,\/d<j)s,i is the eigenvalue of the 
linearized negative drift field at S whose corresponding eigenvector points along <j> S)1 , the direction 
locally parallel to the MPEP. Similarly, Aj_(Z7) = —duup/d&up is the eigenvalue of the linearized 
negative drift field at U whose corresponding eigenvector points along (f>u,2, the direction locally 
perpendicular to the MPEP, and so forth. \\\(U) is the eigenvalue of the linearized negative drift 
field at U that corresponds to the unstable, or 'downhill' direction; it is negative. Eq. (|2TT) . which 
is called the Eyring formula, is a version of the Kramers rate formula that incorporates transverse 
fluctuations. 
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In a spatially extended equilibrium system, where the evolution equation is a stochastic partial 
differential equation (cf. ©), Eq. (OTT) generalizes in a formally (though not computationally) 
straightforward way [ 8 1 . Because of the field-theoretic nature of the model, the linearized dynamics 
at each fixed point have a countable infinity of eigenvalues. The generalization of the prefactor r 
of (ED is 

|A O (0j| I ir=oAn(0.) 

2vr \lir =o |An(0 u )r 1 ; 

where <p s is the stable configuration and U is the dominant saddle. Here A O (0 U ), which is the 
only negative eigenvalue in the sets {A n (0 u )}£L o and {A„(0 S )}£1 O , is the negative eigenvalue of 
the linearization of the negative drift field 5H/5(fi at the transition state <p u . It corresponds to the 
unstable direction. It is easy to see that in the case of a two-dimensional order parameter, in which 
there are only two eigenvalues at each stationary state, (l22ti reduces to the Eyring formula (|2"TT) . 
\\\{U) corresponds to A (0 U ). 

Because both the numerator and denominator of the quotient inside the square root are products 
of an infinite number of eigenvalues (typically with magnitude much greater than unity, as will 
shortly be seen), they may diverge. However, their ratio, defined in a limiting sense, will generally 
be finite. 

It is often not possible to compute in closed form the eigenvalue spectrum of the linearized 
zero-noise dynamics at the relevant stationary points (although an exception will be seen below). 
To employ standard techniques for prefactor computation, the above formula may optionally be 
recast as a determinant quotient of the linearized time-evolution operators at the fixed points ID El 
ED EH- Consider a small perturbation r] about the stable state, i.e., = <p s + r]. Then to leading 
order i] = — A[0 S ]^, where A[0 S ], which is the Hessian operator 5 2 H/5(f) 2 evaluated at = S , 
specifies the linearized zero-noise dynamics at S . Similarly, A[0 U ] specifies the linearized zero- 
noise dynamics at <p u . For notational convenience, we shall write A* and A^ for X n ((f) s ) and A n (0 u ). 
As above, let Aq be the only negative eigenvalue of A[0J, corresponding to the direction along 
which the MPEP approaches the transition state. Then 



r 



l A ol 



detA 



2tt \ |detA[0 u ] 



(23) 



is an alternative expression for T . The two determinants are typically not defined, but their quo- 
tient can be made sense of in any of several ways: as a limit of determinant quotients arising from 
finite-dimensional truncations, for instance. 

In the stochastic Ginzburg-Landau model of ©, it is straightforward to compute the differen- 
tial operators A[0 S ] and A[0 U ]. Linearizing the zero-noise evolution = —SH /5(f) at any stationary 
state = 0o yields 

f] = -A[0 O ] v = ~ {~d 2 /dx 2 + (-1 + 302)] V . (24) 

With the formulas (l23l and (1241 in hand, we shall now work out the Kramers rate prefactor in the 
Neumann case, as a illustrative calculation. 
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5.1 The Neumann- Case Rate Pref actor When L < L c 



If L < = 7r, both the stable and transition states are spatially uniform: <p s = ±1 and <\> u = 0. 
This greatly simplifies the computation of the associated eigenvalues. Note that the eigenvalue 
spectrum is the same at both stable states, by symmetry under <\> i— > — (ft; and because only the 
square (f>l appears on the right-hand side of (l2"3b . 

Linearizing around either stable state yields the operator 

A[(j) s ] = -d 2 /dx 2 + 2, (25) 

and similarly 

k[(j) u \ = -d 2 /dx 2 - 1 . (26) 
The eigenvalue spectrum of A[0 S ], equipped with Neumann boundary conditions, is 

2 2 

A^ = 2 + ^, n = 0,l,2.... (27) 

The eigenvalue spectrum of A[(p u ] is similarly 

7v 2 n 2 

An = -l + -^r, n=0,l,2.... (28) 

As expected, all eigenvalues of A[0 S ] are positive, while the A[<f) u ] operator has a single negative 
eigenvalue (with value —1). The corresponding eigenfunction, which is spatially uniform, is the 
direction in configuration space along which the MPEP approaches 4> u . 

Putting everything together, we find the Neumann-case rate prefactor when L < = n to be 



Z7T 



2 ri 2 > 



noo / Q I 7T z n 
n =oV A + T 1 



2 <n 2 ' 



L 2 



= J_ >^(^) (29) 

23/47rV sinL ' K ] 

where the latter expression follows from the well-known infinite product representation for the 
sine. This diverges as L — > 7r~, i.e., as L — > (L^)~ . In this limit, T ~ const x (L^f — L) -1 / 2 . 
As mentioned in the introduction, this divergence has a simple physical interpretation: it arises 
from the MPEP becoming transversally unstable as L — ► (L^)~ . In the endpoint linearization 
framework, the divergence is caused by the eigenvalue A" tending to zero. The appearance of a 
zero eigenvalue signals the appearance of a transverse soft mode. 



5.2 The Neumann- Case Rate Prefactor When L > L c 

When L > L^, there are two transition states: the nonuniform droplet pair configurations ±0„, 
one of which is shown in Fig. |2j Closed-form computation of the eigenvalues of the linearized 
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dynamics at a nonuniform stationary state is generally not possible. However, there are well- 
known techniques for computing the determinant quotient in ( 123*1) . The transition state (fi u is given 
by ( [""31 . with m = m(L) (for ease of notation, all subscripts will be dropped) determined implicitly 
by (fl4l) . The associated linearized evolution operator, computed from (|24"I) . is 



AW = ~ - 1 + -^r sn 2 \-^= + K(m) 
dx 2 m + 1 V + 1 



m . (30) 



Evaluation of T requires the calculation of the eigenvalue spectrum of this unusual Schrodinger 
operator, or at least the calculation of an associated determinant quotient. It is worth mentioning 
that a similar Schrodinger operator arises in the analyses of Langer and Callan-Coleman. How- 
ever, since they considered only the limit L — > oo, which corresponds torn — ► \~, the 'potential 
energy' in their Schrodinger operators involves hyperbolic trigonometric functions rather than el- 
liptic functions. For a careful review, see Ref. ETl . 

Calculation of the associated determinant quotient is facilitated by the following fact, which 
was first noticed by Gel'fand (c. 1960). Let rj u ^ be a nonzero solution on [0, L] of the homogeneous 
equation 

A[0 u ]t7 = O, (31) 

and let 77^* be chosen similarly. Then if 77^*, 77^* satisfy a Neumann boundary condition at x = 0, 
i.e., 

v' u ,*\o = v' s ,*\o = 0> (32) 

it can be shown that ESIIHI 

det A[0j = v' s JlVu,*\o 

det k[(j) u ] Vu,*\lVs,*\o' 

This reduces the calculation to a series of manipulations of solutions of homogeneous Schrodinger 
equations. 

Consider first the ratio 77^ J l/ Vu,*\o- A solution of the homogeneous equation <f3Tb that satisfies 
the Neumann boundary conditions can be found (see Ref. 1121 ) by differentiating the periodic 
instanton solution (fl3l with respect to m; the result is 



1 + m / 1 \ / m 1 3, . . , nA . 

r]uA z i m ) = / = sn(z\m) — J — ■ r- sn z m (34) 

V 2(1 + m) 1 — m 



[1 — m)^j2m(l + m) 3 

H ; cn(z\m) dn(z\m) 

2m(l + m) 3 



K(m) + (E(m) - E( z\m 

1 — m ^ 



where z = x/^m + 1 + K(m), and cn(- | m) and dn(- | m) are Jacobi elliptic functions [20 1 . 
This yields 

Vu,*\l/vu,*\o = 2 [(1 - m)K(m) - (1 + m)E(m)] . (35) 
The other ratio in (l33l) is easily computed by choosing 77^* = cosh(v / 2x), which yields 

vL\l/vsAo = \/2sinh{V2L) . (36) 
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This is consistent with the numerator of (1291). which was obtained through direct computation of 
the eigenvalue spectrum. Substituting (1331) and (13*61) into 03l yields a formula for the determinant 
quotient. 

Finally, we compute the single negative eigenvalue Aq of A[</> u ], which is associated with down- 
hill motion away from the transition state <j> u . It is easy to check that the corresponding eigen- 
function is 

r] ufi (x;m) = sn 2 ( -f^=f + K(m) m ) - 1 (37) 

\ v +1 / 1 + m — v rrr — m + 1 

In physical terms, this eigenfunction specifies the way in which a droplet pair configuration in 
the Neumann model begins to collapse, due to the boundary between the positive and negative 
droplets, which is a kink initially localized at x = L/2 (see Fig. O tending to move toward x = 
or x = L, where it will be annihilated. The negative eigenvalue itself is 



A" = l Vm 2 - m + 1 , (38) 

1 + m 

which approaches - lasm^ + , i.e., as L — > (L^) + , in agreement with the single negative eigen- 
value of (l28l) . As m — > 1~ (i.e., as L — > oo), Aq — > as const x exp(— Ly/2). This implies that the 
zero-noise motion of the kink domain wall between the positive and negative droplets, which, since 
it leads to annihilation of the kink, is a deterministic 'coarsening' phenomenon (cf. Ref. [ 26 1), pro- 
ceeds only very slowly in the limit of large L. In the same limit, the numerator of the determinant 
quotient diverges only as exp(L/v / 2), so the Kramers rate prefactor tends to zero. 

As an interesting aside, we note that the eigenvalue equation for the operator A[0 U ] of (OUb is 
the spin-2 Lame equation (|2"71. which is a Schrodinger equation with an elliptic potential. This 
potential is periodic, with lattice constant 2K(m). So the spectrum of A[0 U ] has a band structure. 
This observation facilitates the calculation of the eigenvalues of A[0 U ], which can be viewed as 
band edges; and, though we do not supply details here, the calculation of associated determinant 
quotients. The fact that the Lame equation is the eigenvalue equation that governs the stability of 
periodic instanton configurations has previously been noticed by others ll28ll29l . though mostly in 
a quantum context. 

Putting everything together, we find the Neumann-case rate prefactor when L > = n to be 



1 V m? — m + 1 

1 + m 



smh(V2L) (39) 



\ y/2 |(1 - m)K(m) - (1 + m)E(m) 



Compared to the formula (l23l for r , this includes an extra factor of 2, since there are two transition 
states: <p u and — <f> u . As m — > (i.e., as L — > (L^) + ), T diverges as const x (L—L^)~ 1 / 2 . It should 
be noted that as L — > (L*f ) + and L — > (L^)~, the prefactor diverges according to the same inverse 
power law. 

The divergence at L = tt is certainly striking, and leads one to wonder how, when L = n 
exactly, the Kramers rate formula T ~ T exp(— AW/ e) should be modified. Presumably, it should 
be replaced by 

T ~ const x e' a exp(-AW/e) , e -> , (40) 



14 



for some a > 0. The computation of the exponent a and the factor 'const', and of course the 
divergence of the rate prefactor at L = L c for the other choices of boundary condition, will be 
considered elsewhere. 
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